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Plane Couette flow transitions to turbulence for Re ~ 325 even though the laminar 
solution with a linear proflle is linearly stable for all Re (Reynolds number). One starting 
point for understanding this subcritical transition is the existence of invariant sets in the 
state space of the Navier Stokes equation, such as upper and lower branch equilibria 
and periodic and relative periodic solutions, that are quite distinct from the laminar 
solution. This article reports several heteroclinic connections between such objects and 
briefly describes a numerical method for locating heteroclinic connections. Computing 
such connections is essential for understanding the global dynamics of spatially localized 
structures that occur in transitional plane Couette flow. We show that the nature of 
streaks and streamwise rolls can change signiflcantly along a heteroclinic connection. 



1. Introduction 

In plane Couette flow, the fluid between two parallel walls of flxed separation is driven 
by the motion of the walls in opposite directions. Even tho ugh the laminar solution is 



linearly stable for all Re (Reynolds number) as shown by Kr eiss et al. (|l994l ). turbulent 



spot s evolve into large turbulent patches for Re exceeding the modest value of about 
325 ( Bottin et a/.lll998l ). These turbulent patches are sustained by the flow for very long, 
and possibly inflnite, time intervals. From a dynamical point of view, the evolution of the 
velocity fleld corresponds to a trajectory in state space, and indeflnitely sustained motion 
should correspond to invariant sets. Invariant sets in state space have the property that 
a trajectory that starts exactly on such a set stays on that set forever, and a trajectory 
that starts outside that set cannot land on it within a flnite time interval although it can 
approach the invariant set rapidly. Thus a reasonable starting point for understanding 
when and why turbulence becomes sustained in plane Couette flow and other shear flows 
is not the loss of linear stability of laminar flow, which never happens in plane Couette 
flow, but the existence of invariant sets. Equilibria, traveling waves, periodic solutions, 
and relative periodic solutions are all invariant sets. The union of such sets can form 
chaotic saddles or chaotic at tractors, invariant sets which may explain a good deal of 
the dynamics of shear flows ( Schmiegel fc Eckhardtl[l997 ). Thus the numerical compu- 
tation of equi l ibria, traveling waves, periodic solut i ons, and relative periodic solutions 
(lNagatalll99d Il997l: IWaleffeil2QQ3l : fviswanathl [2QQ7I : ICibson et al\\2008^ iHalcrow et al 
20081 Gibson et a/. 12008 ah is a step towards understanding the dynamics of plane Couette 
flow in the transitional regime. 

We use a computational box of extent 0<x<27r/a, — 1<^<1, and < 2: < 27r/7, 



with a = 1.14 and 7 = 2.5 ( Waleffel[2003[ ). where x, z are the streamwise, wall- normal, 
and spanwise coordinates, respectively. Likewise, u^v^w are the three components of the 
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Table 1. Basic statistics for equilibria at Re = 400. EQq is the laminar solution of plane Couette 
flow. The rate of energy input / and the rate of dissipation D are both normalized to be 1 for 
the laminar state. So is the total kinetic energy denoted by E. Eroii/E is the fraction of the 
total kinetic energy in the rolls. The dimension of the unstable manifold is d(W'^), while d{Ws) 
is the dimension of the intersection of the unstable manifold with the 5'-invariant subspace. Ao 
is the eigenvalue with the greatest real part, and Rct is the width of the channel in wall units. 



velocity field. The boundary condition is periodic alon^ x and and no-slip at the walls. 
For comparison, the experimental setup of lBottin et al\ (|l998l ) is about a meter long with 
a separation between the walls of only 7 mm. At the moment, small computational boxes 
are needed to keep the cost of computing invariant sets manageable. Nevertheless, small 
computational boxes are capable of picking up significant aspects of turbulent boundary 
layers and transitional dynamics, perhaps because some of the features of those regimes 
are localized in space. For instance, periodic and relative periodic orbits computed in 
a small box repr oduce the formation and break-up of streaks in the near-wall region 
( Viswanathl 2007 ). Indeed, such solutions show that the spanwise advection of coherent 
structures could be a significant source of the spanwise variation of the root mean square 
value of the streamwise velocity. In addition, turbulent spots, which appear to play a 
key role in transition, are localized in space, although certain aspects of their dynamics 
depend upon the interaction with the external flow (Schumacher & Eckhardt 2001). 

In this article, we report three heteroclinic connections between equilibrium (or steady) 
solutions of plane Couette flow at Re = 400, where the Re is based on half the differ- 
ence in velocity between the moving walls, half the distance between the walls, and the 
kinematic viscosity of the fluid. Basic data for six such solutions, with the first one being 
the laminar solution, is gi ven by Table [T] (for detailed data sets the reader can consult 
Channelflow.org and Hal crowl (|2QQ8l )). The equations of plane Couette flow are un- 
changed by the shift-reflect and shift-rotate transformations defined in Section 3. All 
the equilibria lie in the ^-invariant subspace, which is the space of velocity fields invari- 
ant under both tr ansforniations. The equilib ria EQ^ and EQ2 are called lower and upper 
branch solutions ( NagatalllQQOl : IWaleffell2QQ3[ ) . The heteroclinic connections reported here 
are from EQ3, EQ4, and EQ5 to EQ^; and from EQ^^ to EQ2. 

In the presence of continuous rotation symmetry and discrete reflection symmetry, the 
existence of heteroclinic cycle s follows from the norm a l form of certain codimension-2 
bifurcations ( Kuznetsov 1998 ). Abshagen et al. ( 2QQ4I , 2005 ) have shown that Taylor- 
Couette flow with a stationary outer cylinder undergoes a codimension-2 bifurcation, the 
normal form of which implies the existence of a heteroclinic cycle. That the basic laminar 
solution of this Taylor- Couette flow undergoes a sequence of supercritical bifurcations, 
making it possible to track bifurcations while computing only linearly stable solutions, 
while the transition in plane Couet te flow is subcri t ical is just one difference from our 
work. Notably, the computations of I Abshagen et a/] ( 2QQ4I ) use a domain and boundary 
conditions that match their experimental setup. We do not compute codimension-2 bifur- 
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cations, although we return to that point and the influential thesis of ISchmiegell ( 1999[ ) 
in Section 4. In addition, our computations of heteroclinic connections are explicit and 
make use of the eigenvalues and eigenvectors of the linearizations around the equilibria. 

Instead, our computations rely on the simple principle that an object of dimension k is 
likely to intersect in a stable way an object whose codimension in state space is less than 
or equal to k. At the bottom, this is nothing more than the fact that two submanifolds 
in general position can intersect if the sum of their dimensions is greater than or equal to 
the dimension of the state space (whether they actually intersect is a s ubtle questio n that 
is central to the "structural stability" of er ^odic dynamical systems ( Smald 1967)). For 



an illu stration in the nonlinear setting, see Abr aham fc Shawl (|l992l ). iKevrekidis et al 



( 199Q| ) (see Section 5 of their paper) make elegant use of this principle and of invariant 
subspaces implied by discrete symmetries of the underlying PDE to numerically deduce 
the existence of a heteroclinic connection in the Kuramoto-Sivashinsky equation. Indeed, 
they comment that their work may have implications for shear flows. With regard to the 
heteroclinic connections presented here, it is significant to note from Table [T] that the 
codimension of the stable manifold in the ^'-invariant space (which is equal to d{W^)) 
of EQ^ is less than the value of d{Wg) for EQ^ with i = 3, 4, 5. Thus it is not surprising 
that the unstable manifolds of EQ^ with z = 3, 4, 5 intersect the stable manifold of EQ2 
in a stable way (i.e., robustly with respect to small changes of system parameters). 

All the equilibria in Table [H except EQ5, have well-formed streaks, which means that 
the streamwise velocity has pronounced variation in the spanwise direction. The streaks 
are acco mpanied by streamwise rolls which is the typical situation for boundary layers 
( Kim et a/. 1971). Streaks and streamwise rolls are also found near the ed^es of tu rbulent 
spots (|Dauchot fc DaviaudI Il995l : iTillmarkI Il995l : ISchumacher fc Eckhardtl l2QQlh . They 
could be relevant to the wavelike manner in which the turbulent spots spread to form 
patches. In this regard, we note that heteroclinic connections are important to obtaining 
a global picture of the dynamics in state space. They can be useful for the physical space 
picture as well, as shown by the dramatic change in the balance between rolls and streaks 
along the heteroclinic connection from EQ^ to EQ^. In Section 3, we present a state space 
plot in the manner of lGibson et all (|2QQ86l ) to show how the heteroclinic connections at 
Re = 400 are related to one another. 

2. Finding and verifying heteroclinic connections 

The discretization of the computational box used 32 Fourier points in the x direction, 
35 Chebyshev points in the y direction, and 32 Fourier points in the z direction. Di- 
rect numerica l simulation of plane Couette flow was performed using Channelflow.org 
(lGibsQnll2007r). The equ ilibria listed in Table [T] were found using GMRES-hookstep itera- 
tions (|Viswanathll2007l ) . If the velocity fields of the equilibria are integrated for a certain 
fixed time, they are nearly unchanged. Yet the evolution of perturbations under such 
an integration can be used along with the Arnoldi iteration to determine all unstable 
eigenvalues a nd eigenvectors, as well a set of the least contracting stable eigenvalues and 
eigenvectors ( Viswan ath '2007^ . Such a computation was used to produce the information 
about the unstable manifolds of the equilibria listed in Table [TJ 

The shift-reflect and shift-rotate transformations of a velocity field are given by 



si[u,v,w]{x,y,z) 
S2[u,v,w]{x,y,z) 



[u,v,-w\{x^ L:,/2,y,-z), 
[-U, -V, w]{-x + -y, z - 



(2.1) 



respectively, where and Lz are the periods of the computational box in the x and 
z directions. If either transformation is applied to a trajectory of plane Couette flow, 
one gets another trajectory of plane Couette flow. The space of velocity fields unchanged 
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Figure 1. Plots of distances from the initial (solid line) and final (dashed line) equilibria 
to the velocity field at varying times along the computed heteroclinic connection, (a), (b), (c) 
correspond to the heteroclinic connections into EQ^ from EQg, EQ4, and EQ^, respectively. 



by b oth s\ and 52 is an invariant subspace called the ^'-invariant space in lGibson et al 
([20086). All the computations in this paper are restricted to this invariant space. The 
norm used over velocity fields of plane Couette flow throughout this paper is the square 
root of the kinetic energy, with a normalization that makes the norm of the laminar 
solution equal to 1. 

In a heteroclinic connection, the velocity field of plane Couette flow varies over a time 
(or t) interval infinite in both senses, approaching equilibria as t ^ —00 and as t ^ 00. 
Those are the initial and final equilibria of the heteroclinic connection. Since it is impos- 
sible to integrate over an infinite time interval, our computed heteroclinic connections 
start out in the linearized neighborhood and close to the initial or "out" equilibrium 
Uout, and end in the linearized neighborhood close to the final or "in" equilibrium Uj^, 
after a finite interval of time. For the heteroclinic connections that go from EQ3 and 
EQ5 to EQi, the initial point on the computed heteroclinic connection is a perturbation 
using the two dimensional eigenspace that corresponds to the complex pair of eigenvalues 
with the greatest real part. It is reasonable to look in that space because all except a 
proper subspace of trajectories that originate near an equilibrium point are tangent to 
the leading eigenspace, which is 1-dimensional if the eigenvalue with the largest real part 
is real and simple and 2-dimensional if a simple complex pair. For the same reason, it 
is perhaps surprising that the heteroclinic connection from EQ4 to EQ^ is found using 
an eigenspace corresponding to a complex pair with real part less than that of one other 
eigenvalue. 

Let e^^\ e^^^ be an orthonormal basis for the two-dimensional eigenspace that corre- 
sponds to a complex eigenvalue pair of the equilibrium Uo^t- The span will be tangent to 
the unstable manifold at the equilibrium. We consider the set of velocity fields of plane 
Couette flow that at the initial time T = lie on a circle of radius r: 

11(0)0 = Uout + r(e'^"^^ cos (j) + e*^^^ sin (j)) . 

For a small and fixed value of r, we search for a point on this circle which evolves to 
make the closest approach to another equilibrium, Uj^. Let 

G(0)=niin||u(T)<^-u,„||, 

where u(T)0 is the velocity field that results from evolving the velocity field u(O)0 for 
time T and where the minimizing value of T is the time of the first local minimum 
greater than a certain threshold. The closest approach is the minimum of G{(j)) over 
< (j) < 2iT. Since G(^) is a function of a single real variable, it can be minimized 
using any one of a number of well-known and effective methods. The computation of 
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Figure 2. (a), (b), (c) correspond to EQ^^^ EQ4, EQ5, while (d), (e), (f) correspond to EQ^^^ EQ3, 
EQ5, respectively (there is very little difference in the plots for EQ3 and EQ4). The quiver plots 
in the top row show the streamwise averaged velocity components v and w in the y — z plane. 
The six contour lines of the streamwise averaged u component are equispaced in (umax, —Umax), 
with Umax being 0.44, 0.27 and 0.45 for EQ^, EQ4 and EQ5, respectively, and with the negative 
lines being dashed. The bottom plots show six contour lines of u in the section y — 0. The 
contour lines are equispaced in {umax, —Umax), with Umax being 0.33, 0.22 and 0.18 for EQ^, 
EQ3 and EQg, respectively, and with the negative lines being dashed, lots would be better, 
maybe a stride of 2 in both x and y.] 



heteroclinic connections sketched above uses a first order asymptotic boundary condition 
at the initial equilibrium. For small systems, it is possible to us e an asymptotic boun dary 
condition at the final equilibrium as well. For an example, see lPemmel et all (|2000[ ). 

For the computed heteroclinic connections from EQ3, EQ4, and EQ5 to EQ2, the chosen 
values of r were 0.0001, 0.0003, and 0.0004, respectively. Figure [T] shows data for the three 
computed heteroclinic connections. In each plot of that figure, the solid line is tiny at the 
beginning but rises exponentially while the dashed line is flat. Therefore, we conclude 
that the initial part of each computed heteroclinic connection is in a region where its 
time evolution is governed by the linearization around its initial equilibrium. Similarly, 
we can conclude that the final part is in a region where the evolution is governed by the 
linearization around the final equilibrium. 

To verify the computed heteroclinic connections using another code, it could be nec- 
essary to use three stages. The computed connection from EQ5 to EQ^^, for instance, 
spends about 75 time units near the initial equilibrium and more than 100 units near the 
final equilibrium, as evident from Figure [H Using the data in Table [H one may easily 
estimate that the loss of precision in those two stages is more than 3 digits. As the equi- 
libria themselves are computed with only about 4 or 5 digits of precision, one has to do 
the verification in segments. S uch a verification of Figure [U which was performed using a 
completely independent code ( Vis wanathl 2007), a nd the appHcability of shadowing the- 
orems about numerical trajectories (Palmer 2000l ) leave little doubt that the computed 
heteroclinic connections are real. 
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Figure 3. (a), (b), (c), (d), (e) are plots of the velocity field at t = 50, t = 75, t = 90,t = 100, 
and t — 150, respectively, of the computed heteroclinic connection from EQ5 to EQ^ of Figured] 
The plots are similar to the ones in the top row of Figure [21 with values of Umax being 0.48, 
0.44, 0.36, 0.46 and 0.53, respectively. 



3. Heteroclinic connections at Re — 400 

The top plots of Figure [2] show the correlation between the rolls and the position of 
the streaks. The equilibria EQ^^, EQ3 and EQ4 each have a single counter-rotating pair 
of rolls. The rolls distort the mean flow, which i ncreases with and thus explain the 
position of the streaks in Figures [2^ and b (jKerswell 2005). EQ5 has four counter-rotating 
pairs. From Figure [2f, we see that the mid-plane flow is not at all streaky for EQ5. 

Figure [3] illustrates the manner in which the rolls change in form along the heteroclinic 
connection from EQ5 to EQ^. From Figure [It, it is evident that for t G [75,125] the 
computed heteroclinic connection does not follow the linearized dynamics around its 
initial or final equilibrium. Figure [3] confirms that the rolls change in form within that 
interval . While the coex istence of rolls and streaks in turbulent boundary layers is well 
known (|Kim et a/.lll97ll ), the sort of coalescence of rolls that is observed in Figure [3] is a 
new type of behavior. 

The significance of the heteroclinic connections is that they give a global picture of the 
dynamics, a picture that cannot can be inferred from equilibria alone. To visualize global 
dynamics, it is essential to depict the equilibria and the heteroclinic connections between 
them in state space. The state space of plane Couette flow is infinite dimensional, and 
in the spatial discretization used for computing the heteroclinic connections, it is more 
than 6 x 10^, which is still much too large. Figure [H u ses a 3-dim ensional projection of 
points in that state space, which was introduced bv Gib son et al\ (20086), to depict the 
equilibria and the known connections between them. 

Before explaining the projections used in Figure [H we discuss why that figure can be 
considered a good visualization of the known heteroclinic connections of plane Couette 
flow at Re = 400. It is typical to use projections to construct low dimensional models 
and these models are considered reasonable if they capture 90% of the energy in the 
underlying flow, for instance. Such models use many more dimensions than just three 
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Figure 4. In the state space plot of the three heterochnic connections at Re = 400 above, EQ3, 
EQ4, and EQ5, along with their images under the half shifts t^, Tz and Txz, are denoted by the 
symbols and respectively. The laminar solution EQq (0) is at the origin. EQ2 and its 
images under the half shifts are denoted by O. The meaning of the axes is explained in the text. 

axes, which is all that can be used in a depiction such as Figure IH In addition, if the 
projected velocity field has 90% of the energy, its normwise relative error can be as high 
as 30%. For these reasons, we do not use the amount of energy retained by the projection 
to judge the quality of depictions such as Figure [H 

Instead, we adopt a more geometric point of view. If a curve in a Hilbert space is 
projected onto a finite dimensional plane, it develops artificial cusps or corners at points 
where the tangent to the curve is orthogonal to the plane of projection. For instance, 
the projection of a smooth c urve in to the normal plane at a fixed point on the 
curve has an artificial cusp ([Widdei] [T961I . chapter 3). In Figure IH we see that the 
heteroclinic connections can be wavy but do not have cusps. It is significant that the 
same projection gives a good depiction of all the three heteroclinic connections into EQ^ 
and the heteroclinic connection from EQ^ to the laminar solution, which is shown as a 
thick line. 

To explain the coordinates ai, a2, as in Figured! we define Tx and as follows: 



where and Lz are the periods of the computational box along x and z. In addition, 
'^xz = ^x^z' For each equilibrium that is in the ^'-invariant subspace (|2.ip . one may apply 
and Txz to get three other equilibria that lie in the S'-invariant subspace. In the 
same manner, one may use each computed heteroclinic connection to get three others. 
Only a single copy of each is shown in Figure [H 

Let U2 be the velocity field of the upper-branch solution EQ2, with the laminar velocity 
field subtracted. If e^*^ are defined by 



Tz[u,v,w]{x,y,z) = [u,v,w]{x,y,z^Lz/2), 




Cl(l -^Tx^Tz -^Txz)<^2 
C2(l ^Tx-Tz- rxz)u2 
C3(l -Tx^Tz- rxz)u2 
C4(l - Tx -Tz^ rxz)u2 
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EQ \ I = D E Eroii/E d{W) d{W^) Ao 



EQo 1 1 -0.010966 

EQi 1.710086 0.722516 0.002526 3 1 0.02524949 

EQ2 2.076045 0.634025 0.006357 4 2 0.0441718 

Table 2. The columns have the same meaning as in Table [1] but with the equilibria computed 

at Re 225. 



with C j being normalizing constants, the e^*^ form an orthonormal set ( Gibson et 



200861 ). For a given velocity field of plane Couette flow, the are obtained by subtracting 
the laminar flow from the velocity field and then taking the inner product with e^*\ where 
i = 1,2,3,4. 

The use of the upper branch equilibrium EQ2 to define e^*^ and may appear arbitrary 
and to an extent it is. Heuristically it is a good choice because the computations of 



Gibson et al\ (|20086l ) show that the dynamics of plane Couette flow, including turbulent 
episodes and trajectories that relaminarize quickly, appear to be trapped between the 
unstable manifolds of EQ2 and its three images obtained by applying Tz and Txz and 
the laminar solution. 

4. A heteroclinic connection at Re = 225 

Table [2] gives data for EQq, EQ^^ and EQ2 at Re = 225. By comparing the dimensions 
of the unstable manifolds and their restrictions to the ^'-invariant space, we can infer 
that both EQ^ and EQ2 undergo bifurcations as Re is increased from 225 to 400. The 
dimension of EQ^'s unstable manifold is just 1. By following that unstable manifold, we 
found a heteroclinic connection to EQ2. 

While the dimension of EQ^'s unstable manifold in the /S-invariant subspace is 1, the 
codimension of EQ2's stable manifold in the same subspace is 2. Based on that con- 
sideration alone a heteroclinic connection seems implausible. However, this heteroclinic 
connection is very likely related to a codimension-2 bifurcation. In such a scenario, the 
dimensions of the unstable manifold of the initial equilibrium and of the stable manifold 
of the final equilibr ium must be compared only within the center manifold. 



. eqi 

Schmiegell (|l999) has sy s temat icall y studied bifurca t ions o f the solutions of plane Cou 



Nagatal (|l99Qh and lClever Bussd (|l997h using a representation with 



ette flow found by 

about 1200 modes. He has found heteroclinic connections where the saddle node bifur- 
cation that gives rise to EQ^ and EQ2 is followed soon after by a pitchfork bifurcation 
as Re is increased. The heteroclinic connection reported above is probably of that type. 

To understand this heteroclinic connection better, it could be useful to think of L^, 
the spanwise size of the computational box, as a parameter. In the parameter space 
with Re and Lz as the axes, the saddle-node bifurcations that give rise to EQ^^ and 
EQ2 will form a curve. There will be another curve that corresponds to the pitchfork or 
the Hopf bifurcation. At the intersection of those curves, we will have a codimension-2 
bifurcation. An advantage of realizing a heteroclinic connection using the normal form of 
a codimension-2 bifurcation is that we will get a heteroclinic cycle, not just a heteroclinic 
connection. 

5. Conclusion 

The unstable but recurrent coherent structures observed in turbulent boundary layers 
and in transitional flows are an aspect of turbulent flows. Invariant sets capture some 
features of these coherent structures and their dynamics. While the notion of coherent 
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structures varies with the means used to identify them, the notion of invariant sets is 
much more precise. Compact but hnearly unstable invariant sets in state space (such as 
equihbria, travehng waves, periodic orbits, partiahy hyperbohc tori) are exact solutions 
of the Navier-Stokes equation which correspond to sustained motions of the fluid. 

As a turbulent flow evolves, every so often we catch a glimpse of a familiar pattern. In 
some instances, turbulent dynamics visualized in state space appears pieced together from 
close visitations of equilibria connected by transient interludes. These turbulent interludes 
themselves reflect close passes to other invariant sets in state space, such as unstable 
periodic orbits. Such an approach to turbulence based on a repertoire of recurrent spatio- 
temporal patter ns, which would be periodic or relative periodic or jpits in state space, 
was proposed bv IChristiansen et al\ |l997 ) as an implementation of (119481) 's view 
that turbule nt flows are e rgi:odi c trajectories in state space. A similar approach has been 
suggested bv iNarasimha (|l989l ). who refers to these patterns as molecules of turbulence. 

The heteroclinic orbits that we present here could be the initial steps in charting an 
atlas of the dynamics of plane Couette flow; close passages to equilibria could be identified 
with nodes of Markov graph to give a coarse form of symbolic dynamics, and then these 
heteroclinic cycles would be directed links connecting nodes of the Markov graph. The 
lower branch equilibrium EQ^^, along with the equilibria which connect back to it, appear 
to form a part of the state space boundary dividing two regions: one laminar the other 
turbulent. Turbulent trajectories appear to be trapped between that bou ndary and the 



unstable manifolds of the upper branch equilibrium EQ2, as illustrated bv lGibson et i 
(l2QQ86h . 

The emergence and disappearance of these heteroclinic connections can also be diag- 
nostic. The disappearance of the EQ^ to EQ2 connection is reminiscent of other global 
bifurcations occurring in simpler dynamical systems. For instance, in the Lor enz sys- 
tem a series of such bifurcations occur as the "Rayleigh" number is increased ( JacksonI 
19891 ). For plane Couette flow, such bifurcations could be useful for marking the onset of 



turbulence. 

Future work in this direction should serve to clarify such points. It is still not entirely 
clear what happens at the global bifurcations involved in the creation and annihilation 
of these heteroclinic connections. Furthermore, lists of equilibria and of the heteroclinic 
connections between them found so far should by no means be considered exhaustive. 
Further investigation of plane Couette flow as well as other geometries will most likely 
turn up other dynamically important invariant sets, and more heteroclinic connections 
between them. 
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